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Abstract 

We report on experience with an investigation of the analytic structure of the 
solution of certain algebraic complex equations. In particular the behavior of their 
series expansions around the origin is discussed. The investigation imposes the need 
for an analysis of the singularities and the Riemann sheets of the solution, in which 
numerical methods are used. 
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1 Introduction 



In various problems in theoretical physics we have encountered problems that are formu- 
lated as follows: consider an entire function F of a complex number y, such that 

F(y) ~ y m as y - , (1) 

with nonnegative integer m (in practice, we have met cases with m = 1 and m = 2). The 
task at hand is then to find information about y = y(x) such that 

F(y(x)) = x m . (2) 

In general, both the form of the series expansion of y(x) around x = and the nature 
of its singularities are of interest. Such questions arise, for instance, in the combinatorial 
problem of determining the number of Feynman diagrams contributing to given scattering 
amplitudes in various quantum field theories [|lj , in the computation of the oscillation period 
of nontrivial one-dimensional potential wells D, in the statistical bootstrap model for hot 
hadronic matter (refs. in || and for a recent discussion see 01]), and in renormalization 
theory connected with the 't Hooft transformation ||. An important and interesting 
example, studied in detail in ||, is the so-called bootstrap equation: 

F B (y) = 2y + l- e y , (3) 

which obviously has m — 1. We shall consider functions F of the more general form 

F(y) = P(y) + eyQ(y) , (4) 

where P and Q are polynomials of finite degree dp > and oIq > 0, respectively, with 
real coefficients. As our working example, taken from ||, we shall consider the function 
Fi defined as 

F l (y) = -2-2y + 2e y , (5) 

for which m = 2. In this paper we shall concentrate on the analysis of the Riemann sheet 
structure of those solutions of these equations that have a series expansion around x — 0. 
To determine the asymptotic behavior of these expansions, the nature of the singularities 
will be analyzed numerically. The results are justified by the fact that, in our calculations, 
only finite computer accuracy is required, as we shall demonstrate. 



2 Identification of the Riemann sheets 

As a first step we identify the various Riemann sheets by their value of y(0): the sheet 
labeled s will have y(0) = Y s for that sheet. Obviously, y(0) = is a solution with 
multiplicity m. In general, there will be dp solutions if Q(y) = 0, and infinitely many if 
Q is non-vanishing. It will be helpful if we can identify the Riemann sheet on which pairs 
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s 



Y s /7l 



1 

3 
5 
7 
9 
11 



( 0.0000, 0.0000 ) 
( 0.6649, 2.3751 ) 
( 0.8480, 4.4178 ) 
( 0.9633, 6.4374 ) 
( 1.0478, 8.4490 ) 
( 1.1145, 10.4567 ) 



Table 1: The first few Riemann sheet solutions for Fi(Y B ) = 0. 



(x,y(x)) lie when x is small but nonzero. This is indeed possible, and we shall illustrate 
it using Let us write y = £ + irj with £ and rj real numbers. We are then looking for 
solutions of + irj) = 0, or 



Inspecting the left-hand side of the last equation, we can immediately see that its zeroes 
are quite nicely distributed. We can usefully enumerate them as Im(K.) = u s , where the 
sheet number s takes only the odd integer values ±1, ±3, ±5, .... For positive s, the zero 
u s is certainly located in the interval where smu s > 0, i.e. (s — l)n < u s < sir, and 
U-s = —u s . We have u\ = U-\ = 0, and for increasing s the zero u s moves upwards in its 
interval, until asymptotically we have u s ~ a s — (\oga s )/a s with a s = (s — l/2)n. In Tab.|l| 
we give the values of Y s for Fx, for the first few values of s. Because the values Y s fall in 
disjoint intervals, for small x we need to know y(x) only to a limited accuracy in order to 
be able to identify its Riemann sheet. The only nontrivial case is that of sheets —1 and 1, 
where it is sufficient to consider the complex arguments: for arg(x) — arg(y) = we are on 
sheet 1, for | arg(x) — arg(y)| = 7r we are on sheet —1. Again, limited computer accuracy 
is acceptable here, and for larger m we simply have m different values of the argument, 
distinguished in an analogous manner. Note that of course the labeling of the sheets is 
rather arbitrary: we have chosen the odd integers in order to emphasize that both sheet 1 
and —1 can be considered the principal Riemann sheet. For the bootstrap equation @ it 
is more natural to label the single principal Riemann sheet with y(0) = as sheet number 
zero. 

3 Series expansion 

We want to compute y(x) as a Taylor series around x = 0: 




(6) 




(7) 



n>0 
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Obviously, «o can be chosen as one of the u s above. On principal sheets, with a = 0, we 
also have immediately that ot\ must be chosen out of the m possibilities with a™ = 1. The 
other coefficients must then be computed (algebraically or numerically) by some recursive 
method, which we shall now discuss. 

It would be straightforward to plug the expansion (|7|) into Eq. (|2|) and equate the powers 
of x on both sides, but notice that, for Q non-vanishing, the number of possible products 
of coefficients grows very rapidly, so that the computer time needed to find the first iV 
coefficients grows exponentially with N. As already mentioned in ||, the better way is to 
differentiate Eq. (|2|) with respect to x so that we obtain the nonlinear differential equation 

y'(x) [P'(y)Q(y) + (Q(y) + Q'{y)){x m - P{y))] = mx m - l Q(y) . (8) 

This equation yields a recursion relation involving products of at most dp + dq + 1 co- 
efficients, so that a truncated power series can be computed in polynomial time. As an 
example, for F\ we find the following differential equation: 

y'{x)(x 2 + 2y(x)) = 2x , (9) 

and the following recursion relation: 

olqoli = , 2ao«2 + ol\ — 1 = , 

71— 1 

na a n + (n - 2)a n _ 2 + 2 ^pa p a n _ p = , n > 3 . (10) 

P =i 

We see immediately that y(x) is necessarily even in x if «o 7^ 0, i.e. on the non-principal 
Riemann sheets. In that case, we also see that if a n , n = 0, 2, . . . is a solution, then also 
a* , n = 0, 2, . . . is a solution, where the asterix stands for complex conjugation. This is 
a result of the fact that if y(x) is a solution of Eq. (|5|), then also y*(x*) is a solution. In 
practice, these solutions give the function values on the different Riemann sheets of one 
solution. The analysis of the previous section proves that y s (0) = y_ s (0)* so that the 
solutions satisfy y*(x) = y- s (x*) and the expansion coefficients satisfy 

= (* { n S) Y ■ (11) 

On the principal Riemann sheets we have ao = and a\ = 1 as mentioned, and the two 
solutions on sheet 1 and sheet —1 are related by y~i(x) = yi(—x). For yi(x) we find, 
finally: 



a, 



2(n + l) 



[n 



l)a n _i + 2 ^ POi p a n+ i- 



:i2) 



for n > 2. Using this relation we have been able to compute many thousands of terms. 
The recursion appears to be stable in the forward direction, but we have not tried to prove 
this or examine the stability in the general case. 

In series expansions it is of course always important to know the convergence properties 
or, equivalently, the asymptotic behavior of a n as n becomes very large. In the next section, 
we therefore turn to the singularity structure of y{x). 
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4 Singularities and branches 



In order to find information about the singularity structure of y(x), we employ the tech- 
niques developed in which we recapitulate here. Singularities are situated at those 
values y k of y where 

F'(y k ) = . (13) 

Since F is entire we also know that these singular points must form an enumerable set, i.e. 
we can find, and label, them as distinct points. We shall assume that these singularities 
are square-root branch points, for which it is necessary that 

F"(y k ) + , (14) 

If F" vanishes at y k but F'" does not, we have a cube- root branch point, and so on. If, 
for non-vanishing Q, all derivatives vanish (as for instance when F(y) = e y ) we have, of 
course, a logarithmic branch point. We know that y = — oo corresponds to a logarithmic 
branch point, and it is to remove this to infinity in the x plane that we have required 
dp > 0. In our examples all the singularities at finite x will be square-root branch points. 
The position of the singularity in the x plane, x k , is of course given by 

F(y k ) = x% , (15) 

so that there are m different possible positions, lying equally spaced on a circle around the 
origin. We shall denote them by x k , P with p — 1, 2, . . . , m. Note that, in first instance, it 
is not clear at all whether x k ^ p for certain k and p is indeed a singular point on a specific 
Riemann sheet. Later on, we shall describe how to determine this numerically. For values 
of x close to an observed singular point x kjP we may expand the left-hand and right-hand 
side of Eq. (0) to obtain 

\( y -y h )*F"{y k )~mF(y u )(—-l) , (16) 



2 V 

where we have dropped the higher derivative terms. Very close to the branch point we 
may therefore approximate y(x) by 



Note that there are only two possible values for (3 k ^ p , and each singular point x ktP goes with 
one or the other of these. Again numerical methods will help in determining which one of 
the two is the correct choice. 

We are now in a position to compute the asymptotic behavior of the coefficients a n . To 
find it, we first determine, for a given Riemann sheet, which are the x k)P that lie closest to 
the origin: this gives us the radius of convergence of the expansion of y(x) in that Riemann 
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sheet. We then have to determine those p for which Xk, P is actually a singular point. We 
shall do this numerically, in the way described in the following section. Let us denote the 
set of values of p for which this is the case by V. Now, we may use the fact that 

/ „ (2n-2)\ 

VI -X = 1 - 2_^1 n X , 7 n - 



22n-i( n _ i)! n ! 

n>l v > 

7n ~ -L n-*l 2 + 0(n- 5 / 2 ) , n^oo , 

V47T 



where we have chosen that square root that is real and positive for 1 — x real and positive. 
The asymptotic behavior of a n as n — > oo must therefore be given by 



. (19) 



•1 \ ^ Pk 

/ -> T n 
p£ -p ^k,p 

Amongst other things, this provides a powerful numerical check on the accuracy of the a n 
as computed by the recursive technique. We shall now discuss how the singularity structure 
of our problem can be investigated numerically. 



5 Computer searches for sheet structures 

The main tool we use for our computer studies is a method for taking small steps over 
a Riemann sheet, that is, given the fact that for some value x\ the point y\ = y(xi) is 
determined to belong to a certain Riemann sheet, we perform a small step Ax to a point 
x 2 and find the point y 2 = y(x 2 ) on the same Riemann sheet. Our method to do this is 
nothing but Newton-Raphson iteration: we simply iterate the mapping 

F(y) -xf 

y^y J, * , 20 

until satisfactory convergence is obtained. The starting value for this iteration is just the 
point y\. A few remarks are in order here. In the first place, it must be noted that for this 
method to work, y\ must be in the basin of attraction of y 2 - Since, except at the branch 
points, which we shall expressly avoid, y(x) is a continuous and different iable function 
of x, this can always be arranged by taking Ax small enough. In the second place, the 
accuracy with which y% is actually a solution of Eq. ([|) is not important as long as it is 
in the basin of attraction of y 2 - therefore, there is no buildup of numerical errors in this 
method if we restrict ourselves to just keeping track of which Riemann sheet we are on. 
Finally, problems could arise if two Riemann sheet values of y for the same x are very 
close. But, since F is an entire function, we know that the solutions of Eq. (0) must either 
completely coincide or be separated by a finite distance, any inadvertent jump from one 
sheet to another can be detected and cured by, again, taking a small enough Ax. 

We have applied the following method for detecting and characterizing the various 
singular points. We start on a Riemann sheet s\ at a value x close to zero, and determine 
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y(x) on that Riemann sheet. We then let the parameter x follow a prescribed contour 
that circles a selected would-be singularity Xk >p once (and no other singularities), and then 
returns to the starting point close to the origin. We then determine to which Riemann 
sheet the resulting y belongs. In this way we can find whether Xk, p is, in fact, a singular 
point for the starting sheet, and, if so, which two sheets are connected there. It is also 
possible, of course, to certify the square-root branch point nature of a singular point by 
circling twice around it, and checking that one returns to the original Riemann sheet. 

One important remark is in order here. In our tracking over the Riemann sheet, it 
is necessary that we do not cross branch cuts (except of course the one connected to the 
putative singularity). Since these branch cuts can be moved around in the complex x 
plane, the contour chosen defines the (relative) position of the branch cuts. The sheets 
that are said to be connected at a particular branch cut are therefore also determined by 
the choice of contour. Of course, choosing a different contour will change the whole system 
of interconnected sheets in a consistent manner, so that in fact, given one choice of contour 
and its system of sheets, we can work out what system of sheets will correspond to another 
choice of contour. We shall illustrate this in the following. 

Suppose, now, that Xk jP is one of the singular points on a certain sheet that is closest 
to the origin. We can then follow, on that sheet, a straight line running from x\ close to 
the origin to a point £2 for which X2/xk, P is real and just a bit smaller than one. Since Xk, P 
is by assumption closest to the origin, there is then no ambiguity involved in determining 
which one of the two possible complex arguments of f3k, P we have to take. Thus, we can 
find all the information needed to compute the asymptotic behavior of a n on that sheet. 

6 An example 

Having established the necessary machinery, we shall now discuss a concrete example of our 
method. For this, we have taken the function F\ of Eq. (|5]). It is, in fact, closely related to 
the very well-understood bootstrap equation (|3]): by substituting, in Eq. @, y — > \og2 + y 
and x — > 2 log 2 — 1 — x 2 , we obtain Eq. (|5|). Its Riemann sheet structure, however, is quite 
different, as we shall see. Note that under this transformation, the origin x = 0, y = for 
Fi corresponds to the first singularity in Fg. 

6.1 The singularities 

The values of y(0) on the different Riemann sheets for Fi, namely Y s for s = ±1, ±3, . . . 
have already been discussed above. The singular values yk are simply given by 





(21) 



so that the possible singular points Xk, p satisfy 



x\ p = —Aink . 
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Note that k = does not correspond to a singular point. The positions of the possible 
singularities in the complex x plane are therefore as follows. For positive integer k: 



%k,l — %z k > %h,2 — — lZ k , 
X-k,l = Zk , X-k,2 — ~ z k , 

z k = (l + i)\/2^k . (23) 
At all these various possible singularities, we have 

Pl, P = 8nrk , (24) 
and therefore we may write 

for k > : /3 fe>p = e k>p (l + i)y/Air\k\ , 

for k < : (3 k , p = e k>p {l - i)^\k\ , (25) 

where the only number to be determined is e kjP G {— 1, 1}. It must be kept in mind that 
the value of e depends of course on the sheet: we take the convention that we work on the 
sheet with the lowest number (in absolute value). When viewed from the other sheet, the 
value of e is simply opposite. 



6.2 The Riemann sheet structure 

We now have to discuss how the branch cuts should run in the complex x plane. There are 
two simple options (and an infinity of more complicated ones): in the first option (I), we 
choose to let the branch cuts extend away from the origin parallel to the real axis. This 
corresponds to tracking a contour that, say, first moves in the imaginary direction, and 
then in the real direction, to arrive close to the chosen singularity. The other option (II) is 
to take the cuts parallel to the imaginary axis, so that a contour that does not cross branch 
cuts en route first goes in the real direction, and then in the imaginary direction. Note that 
these two alternatives do, indeed, correspond to different implied relative positionings of 
the branch cuts. In Fig. |l[b we show the contour used in examining singularity a^i under 
option I. The contour starts on sheet number 1 close to the origin (so that y is close to Yi), 
moves upwards and then to the left, circles the singularity once anti-clockwise, and returns 
to its starting point by the same route in order to enable us to determine the resulting 
Riemann sheet number. Fig.|l|.c shows the corresponding path in the y plane. It ends again 
close to Yi so that, for this choice of contour and its induced branch structure (indicated 
in the figure), sheet 1 does not have a branch point at X2,i- Fig. |T[d shows what happens 
if, instead of sheet number 1, we start at sheet number 3: the y track starts then at close 
to Y 3 , but ends up close to Y 5 , so that we conclude that sheets 3 and 5 are connected at 
X2,i- If w e run through the whole contour twice, we get the y track presented in Fig.[l].e, 
where the y track ends up again at Y 3 as expected for a square root branch cut. Under 
option II, we rather use the contour indicated in Fig.|2].a, which first moves to the left and 
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• (3,1) 
' (2,1) 

' (1,1) 


(-3,1) ' 
(-2,1) ' 

(-1,1) ' 


• (-1,2) 
• (-2,2) 
. (-3,2) 


(1,2) . 

(2,2) . 
(3,2) . 



-5 5 

a: The numbering (k,p) of the singularities. 




-2 5 
c: sheet 1, single loop 





-1 



-5 1 

b: x-plane 
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-4 5 

d: sheet 3, single loop 




-4 5 

e: sheet 3, double loop 



Figure 1: The numbering (k,p) of the singularities, and loops around £2,1 under option I. 



"i — 1 r 

1 1 1 
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-5 1 

a: x-plane 



-10 6 

b: sheet 1, single loop 



0.5 3.5 
c: sheet 3, single loop 



Figure 2: Loops around x 2 1 under option II. 
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k 


Xk,l 




%k,2 


Cfc,2 


X-k,l 




X-k,2 




1 


(1,3) 


-1 


(-1,3) 


-1 


(-1,-3) 


-1 


(1,-3) 


-1 


2 


(3,5) 


-1 


(3,5) 


-1 


(-3,-5) 


-1 


(-3,-5) 


-1 


3 


(5,7) 


-1 


(5,7) 


-1 


(-5,-7) 


-1 


(-5,-7) 


-1 


4 


(7,9) 


-1 


(7,9) 


-1 


(-7,-9) 


-1 


(-7,-9) 


-1 


5 


(9,H) 


-1 


(9,H) 


-1 


(-9,-H) 


-1 


(-9,-H) 


-1 



Table 2: Sheets connected at the first few singularities (option I), and the corresponding 
value for e. 
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Xk,2 


Cfc,2 


X-k,l 


e-fc,i 


X-k,2 


C-fc,2 


1 


(1,3) 


-1 


(-1,3) 


-1 


(-1,-3) 


-1 


(1,-3) 


-1 


2 


(1,5) 


-1 


(-1,5) 


-1 


(-1,-5) 


-1 


(1,-5) 


-1 


3 


(1,7) 


-1 


(-1,7) 


-1 


(-1,-7) 


-1 


(1,-7) 


-1 


4 


(1,9) 


-1 


(-1,9) 


-1 


(-1,-9) 


-1 


(1,-9) 


-1 


5 


(1,H) 


-1 


(-1,11) 


-1 


(-1,-H) 


-1 


(1,-H) 


-1 



Table 3: Sheets connected at the first few singularities (option II), and the corresponding 
value for e. 



then upwards. Fig.|2|.b shows the resulting y path, which does not return to Y\ but rather 
to Y5, indicating that under this choice of contour the sheets labeled 1 and 5 are connected 
at X2,i- Fig. |[c shows that, now, sheet 3 is insensitive to this singularity. 

In this way we have mapped the various singularities around the origin. In Tab. § we 
present the pairs of sheets that are pairwise connected at the first few singularities, under 
option I, and the observed value for e, which turns out to be —1 in all cases. We point 
out that at each singularity only two sheets out of all infinitely many are connected. Note 
the somewhat atypical situation at the lowest-lying singularities X\ y ±\ and £_! ±1 . The 
alternative option II results in Tab. |3[ Note that the higher-lying singularities now show 
a sheet structure similar to the lowest ones. In fact, this is the choice that corresponds 
most directly to the analysis of the sheet structure of the bootstrap equation in [[| , with of 
course the extra complication in the fact that the bootstrap equation (|3]) has m = 1 while 
for Fi, m = 2. Note that, once again, e = —1 in all cases. 

6.3 Asymptotic behavior of the series expansion coefficients 

We shall now illustrate how the information on the Xk )P and /3k )P allows us to compute the 
asymptotic behavior of the series expansion coefficients a n . 
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2 - 
1 ■ 

0- 
-1 ■ 
-2 ■ 




logn 



Figure 3: r n , defined in Eq. (ESp, as function of logn. 



First Riemann sheet. In this sheet, the singularities closest to the origin, and their 
corresponding /3's are 



xi t i = v 47r exp(3i7r/4) , /3i 5 i = v 87rexp(— 3i7r/4) , 

x_ 12 = v / 4Texp(-3i7r/4) , /3_ 1>2 = v / 8^exp(3m/4) . (26) 
Using Eq. (O), we see that the asymptotic form of the coefficients on sheet 1 is given by 



c/ 1 ) « asym 

n n J 

• 1 . » . .-. 1 /2 n ' 



n 3 / 2 (47r)"/ 2 

/- /3n7T 37T \ I u // — i/; i 

-v / 2cos(^ + T ) = <| , (27) 



with integer p. In Fig.|3|. we have plotted the observed behavior of 
/f4vrW 2 n 3 / 2 \ 




on the first Riemann sheet, against logn. The coefficients clearly converge to the computed 
behavior, and we can even distinguish that the leading corrections go as n~ 5 ^ 2 ; the four 
separate lines that emerge are just the four different forms of c n . The series expansion for 
Riemann sheet —1 are simply obtained from 

= R n «i 1} • (29) 
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Higher Riemann sheets. We first consider positive sheet label s = 3, 5, 7, . . . and put 
k = (s — l)/2. We then have 



Xk,i = ~Xk,2 = v47r/cexp(i7r/4) , (3 k ,i = fa, 2 = (1 + i) V^rk . (30) 

As we have already seen a n vanishes for odd n, and for even n we have the following 
asymptotic form: 

(s) 2(l + i)Vk +1 w 2(l-i)Vk 

~ (4**)* (_) ' " 4p+2 ~ ' (31) 

for integer p. For negative s, we use Eq. (pTT|), which also holds asymptotically. 



7 Conclusion 

We have shown how numerical tracking can be helpful in the investigation of the Riemann 
sheet structure of the solution of certain algebraic complex equations. Furthermore, we 
have shown how the series expansions around the origin on the different sheets and the 
asymptotic behavior of their coefficients can be determined. The results of the numerical 
analyses have been justified by the fact that only finite computer accuracy was required in 
the specific calculations. 
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